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Abstract. For a prescribed pair of quantum states \ipx) an d \iI>f) we establish 
an elementary derivation of the optimum Hamiltonian, under constraints on its 
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(N 



| eigenvalues, that generates the unitary transformation \ipi) — > \ipF) in the shortest 

duration. The derivation is geometric in character and does not rely on variational 
1"**- , calculus. 
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Recently Carlini, et al. pQ considered the following problem: What is the optimum 
choice of the Hamiltonian, under a given set of constraints, such that the transformation 
between a designated pair of quantum states is achieved in the shortest possible time? 
Evidently this question is of relevance to the implementation of various algorithms 
in quantum computation (see, e.g., |TJ Ej and references cited therein). Two specific 
examples for the constraints on the Hamiltonian are considered in pQ, and the optimum 
O" 1 solutions are obtained using the method of variational calculus. 

The purpose of this paper is to show that analogous results can be obtained more 
directly by use of the symmetry properties of the quantum state space, hence avoiding 
the use of variational calculus. Our approach is closely related to the idea considered 
in where an elementary derivation is provided for the minimum time required to 
transform one quantum state into another for a given Hamiltonian. The idea here is to 
reverse the argument to find the optimum choice of the Hamiltonian that achieves the 
transformation in the minimum time. 

Consider a Hilbert space 7i n+1 of dimension n+1, and assume that an arbitrary pair 
of initial and final states and \i[>f) are specified. The task is to find the Hamiltonian 
H on 7i n+l that generates the unitary transformation \tpi) — > \ipp) = e lHr ^ h \tpi) in 
shortest possible time r. Clearly, if the differences between the eigenvalues of the 
Hamiltonian are allowed to take large values, then the value of r can be made very 
small. This is because the 'speed' of a unitary evolution is proportional to the energy 
uncertainty (the so-called Anandan-Aharonov relation 0j). As a consequence, if the 
differences between eigenvalues can be made large, the energy uncertainty can also be 
made large. Hence we impose the constraint that the difference of the largest and the 
smallest eigenvalues of H be bounded by a constant. 
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Figure 1. Optimum state transformation. In Hubert space H n+ one wishes to 
unitarily transform the initial state into the final state \4>f) in the shortest possible 
time. In 7Y" +1 there is a unique two-plane Tt 2 that contains the two endpoints of the 
vectors \tpi) and IV'-f), and the origin. In projective terms this plane corresponds to a 
complex projective line in the space V n of pure states, and the two vectors \ipi) and 
\iPf) determine a pair of points on this line. The geodesic curve that joins these two 
points lies on this complex line, which in real terms is just a two-sphere S 2 . Given 
a pair of points ipi and ipF on S 2 there is a unique great circle arc passing through 
these points. The most expedient transformation \ipj) — > \iPf) is thus obtained by the 
rotation of S 2 around the axis that is orthogonal to the hemispherical plane containing 
■0/ and if) p. The axis of rotation, in particular, corresponds to a pair of orthogonal 
states \E±). The Hamiltonian that generates this rotation therefore takes the form 
H = X + \E + ){E + \ + \-\E-){E-\ for a pair of real parameters A±. 



For the analysis of a problem of this kind it is useful to work directly with the 
space of rays through the origin of 7i n+1 . This is just the complex projective space V n 
of dimension n; each ray \ip) e 7i n+1 then corresponds to a point ip G V n . Thus V n can 
be thought of as the space of directions in Ti n+1 . Now given a pair of points if)j, ipF G V n 
corresponding to the states \ipi) and \ipp) in 7i n+l we can join these two points by a 
line. The points on this line correspond to all possible linear superpositions of the states 

and \iPf)- That is, the (complex) line in V n corresponds to the two-dimensional 
subspace of H n+1 spanned by the two vectors \ipi) and \ipf) (see FigureHJ). In real terms 
the complex line in V n corresponds to a two-sphere (the so-called Bloch sphere) S 2 , and 
the two states thus correspond to a pair of points on the surface of this two-sphere. 

It is evident that there is a unique geodesic curve on S 2 that joins ipj and tjjp] this 
is just the great circle arc passing through these two points (cf. [HUH]). Therefore, the 
unitary transformation that takes the state into in the smallest possible time 
is achieved by a rotation of S 2 around the axis such that the geodesic curve joining 
ipi and ipF constitutes the equator associated with that axis. There are infinitely many 
other unitary transformations that achieve the transformation \ipf) — ► \i/)f), however, all 
these transformations will require longer times to be realised because the corresponding 
trajectories are not geodesic curves. 

We thus proceed to determine this axis of rotation. To this end let us write \%pi) for 
the state orthogonal to the initial state that is contained in the two-dimensional 
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span of the initial and final states in 7i n+1 (i.e. the antipodal point on S 2 ). Then the 
final state can be written in the form 

|V>f) = cos \e\^i) + e i( ^ /2) sin §0^). (1) 

Since both |^/) and \ipp) are prespecified, the values of the two parameters 9, <fi are 
known. Our objective now is to find the axis defined by a pair of antipodal points on 
S 2 for which and \^f) lie on the equator (see Figure C}. Since \ipi) and lie on 
the equator associated with the (E + , i?_)-axis, conversely the two states \E+), \EJ) lie 
on the equator associated with the (if)j, ^/)-axis. Hence these states can be expressed 
as equal superpositions of \ipi) and \4>i): 



E +) = 75 (IV*> + e 1 ^/)) and |£L> 



i 

v/2 



Solving (JTJ for and substituting the result into (J2J) we obtain 



and 



'*-> = 71 L 
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sin 
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cos \9 
sin \9 



\1>l) + 



sin 1 6* 



sin i# 



(2) 



(3) 



(4) 



These states thus determine the axis of rotation that we are seeking. 

Now the unitary rotation that gives rise to the rotation of the two sphere about the 
axis \E + ) and \E_) is generated by the Hamiltonian 



H = \ + \E + )(E + \ + \_\E_)(E_\ 

for some choice of real parameters A + ^ A_. Substituting 
express this Hamiltonian in terms of the two input states: 
A + + A_ 



and 



(5) 

into (jSJ) we can 



H 



2 sin 2 \9 



+ 



+ 



A, 



A, 



cos 1 9 



sin \9 



sin 2 \9 



— l- 



sin ^ 



cos |0 
sin 2 p 



A_ 
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A 



1 



sin ~0 



+ 



cos ~0 



sin 2 ±0 



+ 



sin 



cos \9 
sin 2 p 



(6) 



Because the Hamiltonian in standard quantum mechanics is defined up to an overall 
additive constant, without loss of generality we may set A + — A_ = £, and hence, 
A + = — A_ = £/2, for some real parameter £. It then follows at once from (jUJ) that 

1 1 



2 sin \9 



2 sin \9 



(7) 



Finally we shall impose the constraint that the difference of the largest and the smallest 
eigenvalues (here there are only two) of the Hamiltonian be given by 2u. Since the 
eigenvalues of H in (J7J) are ±£/2 sin |0 we have ui = £/2 sin |0. 

More generally, we may consider time dependent Hamiltonians. However, because 
of the constraint on the difference of the eigenvalues, the parameter uj cannot vary in 
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time. As a consequence the only time-dependence that can be introduced here is that 
associated with the overall magnitude of the Hamiltonian, which in itself does not affect 
the dynamics. Letting h(t) denote this gauge term and 1 denote the identity operator, 
the optimum choice for the Hamiltonian can thus be written as 

H = iu#/)(ViH - M<M(^I + h(t)l. (8) 

This is the main result obtained in pQ. We emphasise that this result is obtained here 
from the symmetry properties of the quantum state space, essentially only requiring the 
use of elementary trigonometry. 

As noted above, the time it takes to achieve the transformation —>■ \ipF)i 
under the unitary evolution generated by the Hamiltonian (jHJ), can be determined from 
the Anandan-Aharonov relation [I], which states that the 'speed' of the evolution of 
a given quantum state is given by 2h~ 1 AH, where AH is the standard deviation of 
the Hamiltonian. Note that the energy variance is a constant of motion. Therefore, 
calculating the standard deviation of (jEJ) in the state, say, we deduce that 

AH = u sin -9. (9) 

On the other hand, the separation of the two states and is just the angle 9. 
We thus find that 

r = j^r 9 . (10) 

2uo sin 2" 

Alternatively, the time required for achieving the transformation can be determined 
more explicitly as follows. We take the Hamiltonian (jHJ) and use it to calculate the 
time-dependence of the state explicitly as 

cos (h^ut sin 1$) : — \- sin {fi~ x ujt sin |#) 

sin ^ u 

+ -^ T -sin(7r 1 wtsin|fl) \ip F ), (11) 
sin 2^ 

where \ip{0)) = Evidently the coefficient of \ipi) in the state \i[>(t)) first vanishes 

at time t = h9/2ujsm while at that time the coefficient of becomes unity. 

In the foregoing material we have considered the case for which there is only 
one constraint on the Hamiltonian, namely, that the difference of the largest and the 
smallest eigenvalues be a constant. In a more realistic setup, however, there can be 
further constraints to limit the allowable operations. Although the use of variational 
calculus suggested in [T] is quite effective in general, it should be evident that within a 
given context, the determination of the optimum Hamiltonian that achieves the desired 
transformation simplifies considerably by taking into account the symmetries of the 
relevant state space. 
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